Enhanced spontaneous emission from inhomogeneous ensembles of quantum dots is 

induced by short-range couplings 
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We study theoretically the spontaneous emission from an inhomogeneous ensemble of quantum 
dots in the weak excitation limit. We show that collective, superradiance-like effects lead to an 
enhanced emission rate in the presence of sufficiently strong coupling between the dots in agreement 
with experimental observations, which means that the quantum dot sample cannot be treated as an 
ensemble of individual emitters. We demonstrate also that the collective behavior of the quantum 
dot system relies on short-range interactions, while long-range dipole couplings are too weak to have 
any impact on the emission dynamics for a system with a realistic degree of inhomogeneity. 
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The unique atomic-like optical properties of semicon- 
ductor quantum dots (QDs) allow one to transfer ideas 
and technologies from atomic systems to solid state struc- 
tures and devices. Such quantum optics concepts success- 
fully implemented in QD systems range from laboratory- 
scale coherent optical control experiments [l| to commer- 
cially available devices, like QD lasers Q. While similar- 
ities between atoms and QDs may be very useful, these 
two systems differ essentially with some respects. One 
of these differences is the inhomogeneity of the QD char- 
acteristics due to inevitable randomness of the physical 
properties of these structures. This becomes important 
when collective interaction of QDs with their environ- 
ment is essential, e.g., in lasers or some quantum informa- 
tion devices [3J . Understanding the collective evolution of 
non-identical, coupled emitters can also be beneficial to 
the current study of other semiconductor plasmonic 
[B| , hybrid Q or atomic Q systems where a similar inter- 
play of collectivity, inhomogeneity, and interactions plays 
a crucial role. 

In an experiment performed on a dense ensemble of 
CdSe QDs [8j , the decay of luminescence was shown to 
accelerate as the number of emitting dots increased. This 
means that, in their interaction with the radiation field, 
the QDs in the ensemble are not independent objects. 
The existence of such a cooperative effect for QDs is re- 
markable as the ensemble broadening of the transition en- 
ergies is 3-4 orders of magnitude larger than the radiative 
broadening of the QD emission line. The experimentally 
observed effect was quite strong already for fewer than 
100 QDs so that an accidental spectral overlap between 
two or more dots was clearly very unlikely. The cooper- 
ative effect was therefore attributed to the existence of 
long-range (LR) coupling between the dots. 

Subsequent theoretical analysis 0-11 confirmed that 
coupling between non-identical QD emitters stabilizes 
the collective effects in the spontaneous emission. How- 
ever, in order to overcome the inhomogeneity the cou- 
pling must be at least comparable to the transition en- 
ergy mismatch between the emitters [j| [l(|, while the 



LR dipole (Forster) couplings do not exceed a fraction of 
meV even for dots separated be a few nanometer distance 



12] and drop down to about 1 /ieV for inter-dot separa- 



tions around 30 nm, which is the average value for the 
experimentally studied sample Q. These values are at 
least two orders of magnitude lower than the transition 
energy inhomogeneity which is on the order of a few tens 
of meV. Thus, the LR dipole couplings seem unlikely to 
underlie the enhanced emission. It is known, however, 
that also other couplings of different nature exist in QD 
ensembles plU^ . 

In this paper, we present the results of theoretical mod- 
eling of the spontaneous emission from an inhomogeneous 
ensemble of QDs and clarify the origin of the observed 
enhanced spontaneous emission. We extend the stan- 
dard model of identical emitters flEl [l6| by including 
ensemble inhomogeneity and simulate the evolution of 
a planar ensemble of non-identical, randomly distributed 
QDs in the weak excitation regime. Our analysis con- 
firms that the cooperative interaction of the QDs with the 
surrounding electromagnetic (EM) vacuum can lead to 
an increased spontaneous emission rate. We show, how- 
ever, that the dipole coupling between the QD emitters 
is insufficient for the appearance of cooperative behav- 
ior of these strongly inhomogeneous systems. Enhanced 
emission appears only if one includes short-range (SR) 
couplings between the dots, which may be due to a com- 
bination of tunneling (wave function overlap between the 
neighboring QDs) and Coulomb correlations. 

We model a system of N QDs located at points r a . 
The typical energy distance to p-shell exciton states in 
CdSe/ZnSe dots is at about 100 meV [T^J, much larger 
than the inhomogeneous broadening of the s-shell states, 
which allows us to restrict the discussion to the funda- 
mental transition. In addition, we fix the polarization. 
Hence each of the dots is represented as a two-level sys- 
tem with the basis states |0) Q and |1) Q , where a labels 
the dots. The fundamental transition energy in the dot 
a is E a . The dots are coupled by SR couplings de- 
scribed effectively by the coupling potentials ii~ a J. The 
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Hamiltonian for the EM interactions is transformed to 
the dipole form using the Power-Zienau-Wooley trans- 
formation Il8f and reads 
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where the EM field is represented by the displacement 
operator 
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Here b k \ is the annihilation operator for a photon with 
wave vector k and polarization A, uik is the photon fre- 
quency, e k \ is the unit polarization vector, a a — (|0)(1|) Q 
is the transition operator for the dot a, do is the inter- 
band matrix element of the dipole moment (we assume 
that all the interband dipoles are identical), c is the speed 
of light in vacuum, £o is the vacuum permittivity and e r 
is the dielectric constant of the semiconductor medium. 

The evolution equation for any electronic operator Q 
can be obtained along the lines worked out in Ref. fl9| : 
One writes down the equations of motion for the elec- 
tronic and photonic operators, eliminates the latter, ne- 
glects the off-resonant terms and radiation-induced en- 
ergy shifts, and performs the Markov approximation. It 
is assumed that no external fields are present. The re- 
sulting evolution equation for the average value of the 
electronic operator Q has the form 

(G) = \ E £ « ( [ <J a' J ai Q] ) + £ ^0 { , Q] ) 

+ J2 T ^^ a Qa p -^{aia p ,Q}y (1) 

where Q. a p = 0,^p + fij^ and 

^(r^d) _ TG a p(k r a p), T a p = TF a p(k r a p). (2) 

Here r a p = rp — r a , T = \do\ 2 ko/(3irso£ T fi) is the sponta- 
neous emission (radiative recombination) rate for a sin- 
gle, isolated QD, ko = nE/(Hc) is the average resonant 
wave vector in the dielectric medium with the refractive 
index n, 
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and du and d± are the components of do = do /do par- 
allel and perpendicular to f a p, respectively. The first 



expression in Eq. @ describes the dipole couplings be- 
tween the interband dipoles confined in the QDs (referred 
to as Forster coupling or dispersion forces). We as- 
sume that only heavy hole excitons are involved, so that 
do = d (l,±«,0) T /v/2, hence |J||| 2 = \dj_\ 2 = 1/2 since 
the QDs are distributed in the xy plane. An admixture 
of light hole states would reduce the dipole coupling at 
short distances. However, typical admixture of at most 
several per cent [2(| will only bring a small quantitative 
correction to our results and will not affect the qualita- 
tive conclusions of this paper. 

In order to find the evolution of the total number of 
excitons, Nx = J2a( a a a a) < we use Eq. ([T]) to find the 
equations of motion for the quantities x a p = {cr a crp). 
These quantities are dynamically coupled to higher order 
terms of the form (a^a^cr^as) ■ In general, in an inhomo- 
geneous ensemble, there are no constants of motion in 
the iV-dot dynamics and the system is described by the 
2 2N elements of the general density matrix, which is be- 
yond the simulation capabilities already for several dots. 
Here, we will restrict the discussion to the low excitation 
case, with at most one exciton present in the QD ensem- 
ble (in this case, collectivity results from derealization of 
this single excitation over many emitters jointly interact- 
ing with the field). Then, the normally ordered higher 
order terms vanish and one obtains a closed system of 
equations of motion in the form 

x a p = — (e Q — ep) x a j3 + i ^ (r2 7Q x 7/ 3 — Qp 7 x ai ) 

7 

— — ^ ' (T~fa x ~fp + Ff3~fXa~f) j (3) 



where it is assumed that f2 77 = 0. 

The photon detection rate for a hypothetical ideal de- 
tector collecting from a solid angle df2 around a point R 
at a large distance from the sample is proportional to the 
field correlation function 211, 



J(R)dQ = 



2CH 2 



D^(R,t)-D^(R,t))dQ, (4) 



where are the positive and negative frequency parts 
of the field. The field originates from the interband 
dipoles and can be related to the transition operators 
<7„. In the Markov limit one finds 1 191 
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"• a (ai(t)ap{t)). (5) 



a/3 



For a spectrally resolved detection, one has to replace the 
field operators in Eq. (U) by spectrally filtered operators 
Z)(±)(i) = J dt'f(t - t')D^{t'), where f(t) is the in- 
verse Fourier transform of the filter function /(w). This 
results in an analogous replacement for the a a , ap oper- 
ators in Eq. ([5]). The two-time correlation functions are 
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FIG. 1. The time dependence of the exciton occupation for 
a QD ensemble on a 0.03 fim 2 mesa (30 QDs on the aver- 
age), (a) With LR dipole couplings only. Solid (red) line: 
homogeneous system. Dashed lines: three different inhomo- 
geneous systems with the standard deviations of the funda- 
mental transition energies a as indicated, (b) With additional 
SR couplings (solid red line) and with the dipole coupling ar- 
tificially enhanced by a factor of 400 (dashed green line) for 
a = 28 meV. In both figures, the dotted (gray) line shows the 
exponential decay with the rate F characteristic of a single 
dot. The results are averaged over 20 random QD distribu- 
tions (the same for each case). 



calculated by means of Eq. ([T]) and the quantum regres- 



sion theorem [21| . starting from the equal time average 
known from the solution to Eq. (J3J . For a milli-electron- 
volt filter width, the relevant time differences are on the 
picosecond scale, much shorter than the exciton life time, 
so that the dissipative part of Eq. ((T|) is neglected. 

To model the experimental situation Q, we consider 
ensembles of quantum dots placed on square mesas. The 
spectral properties of the dots are modeled by a Gaussian 
distribution of their transition energies with the mean E 
and variance a 2 . The dots are located at random in the 
sample plane, with the restriction that the inter-dot dis- 
tance cannot be smaller than 10 nm. The mesas are "cut 
out" from a larger sample, so that the number of dots 
fluctuates. At the initial time, the QD ensemble is sup- 
posed to be coherently excited into the delocalized state 
|^(0)) = (l/jV) 1/2 £ Q cri|0), where |0) is the system 
ground state. For each set of parameters, the results are 
averaged over many repetitions with different, randomly 
generated QD distributions. In the simulations presented 
here, we assume 1/T = 390 ps, n = 2.6, E = 2.59 eV and 
the QD surface density v = 10 11 QDs/cm 2 . For the SR 
coupling between the QDs we choose the simplest expo- 
nential model, £l a J = Vq expj— r a p/ro], with the ampli- 
tude Vq — 5 meV and the spatial range ro = 15 nm, 
chosen to fit the experimental data. 

The decay of the exciton population in a system of 
30 QDs on the average with various degrees of inhomo- 
gencity and coupled by interactions of different kinds is 
shown in Fig.[TJ If only LR dipole couplings are included 
[Fig. [2a)] the evolution of a system of identical dots is 
markedly non-exponential with a long stage of strongly 
accelerated emission (red solid line). However, this ef- 
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FIG. 2. (a) Dependence of the average exciton life time on the 
"mesa" size used in the simulation for an ensemble of dipole 
coupled dots (red squares) and SR-coupled dots (green cir- 
cles). The error bars show formally the values of the standard 
deviation of the results obtained for various QD distributions; 
however, the actual distribution of results is strongly asym- 
metric with a tail towards shorter life times and no results 
above the single dot exciton life time of 0.39 ns. (b) Average 
exciton life time as a function of the mesa size for two hy- 
pothetical QD ensembles: with dipole forces scaled up by a 
factor of 400 (blue triangles) and with SR couplings but emit- 
ting to separate reservoirs (brown diamonds) . The results are 
averaged over a number of repetitions ranging from 10 4 for 
the small systems to 100 for the largest system studied. 



feet vanishes very quickly as soon as the inhomogeneity 
of the fundamental transition energies comes into play. 
As can be seen in Fig. Q] (black dash-dotted line), already 
for a = 0.1 meV, the decay almost exactly follows the 
exponential evolution of a single dot, with only a very 
short initial period of enhanced emission. 

The results are completely different if a sufficiently 
strong SR coupling is taken into account [Fig. [IJb), red 
solid line]. Now, the exciton occupation decays with an 
enhanced rate even for the realistic value of a = 28 meV 
0. Moreover, although the time dependence is not 
strictly exponential, the emission rate is much more con- 
stant over the two decades of occupation change shown. 

Although one could expect that the LR dipole interac- 
tion should lead to qualitatively different dynamics than 
the SR coupling this turns out not to be the case for the 
inhomogeneous QD system. As shown in Fig. db) (green 
dashed line), the evolution of the exciton occupation in 
the hypothetical case of dipole couplings magnified by 
a factor of 400 is almost identical to that observed for 
SR couplings for the parameters of these couplings cho- 
sen here. The reason is that in any case the coupling 
is important for the spontaneous emission only when its 
magnitude is at least comparable to the energy mismatch 
between the emitters Q. For a realistic, rather strongly 
inhomogeneous ensemble of QDs, only coupling between 
relatively close QDs is strong enough to have any impact 
on the system dynamics. Hence, the properties of the 
interactions on longer distances are irrelevant. 

As pointed out in Ref. [8[, a fingerprint of collective 
spontaneous emission is the dependence of the exciton 
decay rate on the number of emitters. In Fig. HJa), we 
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show the simulation results for this dependence in the 
case of a QD sample coupled only by the dipole forces 
(red squares) and in the presence of the SR coupling 
(green circles). As expected from the discussion pre- 
sented above, there is no noticeable decrease of the ex- 
citon lifetime in the case of LR EM coupling. However, 
including the SR coupling leads to a considerable acceler- 
ation of the spontaneous emission to a degree comparable 
to that observed in the experiment @ . A similar result is 
obtained in the hypothetical case of artificially magnified 
dipole couplings, as shown by blue triangles in Fig.^b). 

Interestingly, for both kinds of couplings, the effect of 
a variation in the QD surface density is quite similar: We 
have checked that decreasing the density by a factor of 4 
reduces the increase of the decay rate) for a 0.03 /Ltm 2 mesa 
to 23% and 25% of the original value for the SR and en- 
hanced LR coupling models, respectively. Increasing the 
density by the same factor enhances this collective effect 
by similar factors of 3.3% and 2.6%, respectively. At the 
same time, even in a sample with increased density, the 
normal dipole coupling is still not sufficient to induce a 
noticeable enhancement of the spontaneous emission. 

The results of simulations show that strong SR cou- 
pling between neighboring QDs is crucial for the enhance- 
ment of the spontaneous emission, while the LR dipole 
coupling encompassing the whole sample is irrelevant. It 
must be noted, however, that the collective nature of the 
interaction between the QDs and the EM field is still es- 
sential for the observed effect. As follows from Eq. ([1}, 
collective EM coupling affects the system dynamics in 
two ways: It not only mediates the dipole interactions 
described by the coupling constants but also leads 

to the appearance of interference terms T a p for a ^ (3 in 
the dissipative part. These terms are absent in the hypo- 
thetical case of QDs emitting to separate reservoirs (but 
still coupled by interactions) The brown diamonds 
in Fig. [2] show the exciton life times for such an artificial 
system of QDs with additional SR interactions and cou- 
pled to separate reservoirs. Clearly, no enhancement of 
the spontaneous emission is observed. 

Another experimental fingerprint of collective emission 
is the dependence of the luminescence decay time on the 
detection energy Q: Since in the spectral range around 
the maximum of the ensemble photoluminescence (PL) 
more dots contribute to the emission the collective effects 
in this energy sector are stronger and the decay is faster. 
Our simulations, shown in Fig. [31 show that also this 
feature is reproduced in a model with SR coupling. Here, 
we have simulated the measured signal from an ensemble 
with a = 28 meV, using a Gaussian spectral filter /(cj) ~ 
exp{— (l/2)[(w — oj )/Auj} 2 } with HAuj — 3 meV, and hwo 
at the PL maximum or 40 meV away, for two sizes of the 
simulated mesas. As can be seen, in both cases the PL 
decay at the PL maximum is indeed faster than in the 
tail of the inhomogeneous distribution. 

Our simulations of spontaneous emission from ensem- 
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FIG. 3. Time dependence of the luminescence intensity (nor- 
malized to the initial value) for detection at the photolumi- 
nescence maximum (solid red lines) and 40 meV away from 
the maximum (dashed green lines): (a) for a 0.01 /mi 2 mesa; 
(b) for a 0.1 fj,m 2 mesa. The results are averaged over 100 
different QD distributions. 



bles of dipole- and SR-coupled QDs show that the emis- 
sion rate in such a system can indeed be increased due to 
collective coupling of the emitters to the EM field. How- 
ever, in view of relatively large inhomogeneity of the QD 
transition energies, sufficiently strong coupling between 
the dots is needed to stabilize the collective nature of the 
emission. For typical inter-dot separations, fundamental 
dipole interactions are too weak to play an important role 
in the emission kinetics. However, the presence of SR in- 
teractions, which may be due to some kind of electronic 
coupling between the dots in combination with Coulomb 
correlations, leads to enhanced emission in quantitative 
agreement with the experimental results Q . While such 



a coupling is likely to exist in QD samples [13j, its ex- 
act microscopic nature is neither important for the emis- 
sion dynamics nor can be inferred from it. In any case, 
the presence of collective effects in the emission means 
that QDs are not necessarily independent emitters, as 
usually assumed when modeling their optical properties. 
Whether the increase of emission rate at low excitation 
is a signature of the true superradiance that could re- 
sult in a delayed outburst of radiation under strongly 
inverted initial conditions 22J remains an open question. 
Also the evolution of inter-dot coherence in the process of 
carrier trapping and relaxation, which would account for 
the experimentally observed differences between emission 
under resonant and non-resonant conditions, remains an 
interesting problem for further study. Finally, in the 
course of further research it may me interesting to clarify 
whether this kind of cooperative behavior can be respon- 
sible for radiative transfer of excitation between remote 
dots [H. 
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